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Abstract The perturbation of the free rigid rotator by the trigonometric Scarf potential is shown 
to conserve its energy excitation patterns and change only the wave functions towards spherical 
harmonics rescaled by a function of an unspecified parity, or mixtures of such rescaled harmonics 
of equal magnetic quantum numbers and different angular momenta. In effect, no parity can be 
assigned to the states of the rotational bands emerging in this exotic way, and the electric dipole 
operator is allowed to acquire non-vanishing expectation values. 
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1 Introduction 

The energy spectra of molecules consisting of two atoms are as a rule well understood in terms of 
equidistant vibrational excitations characterized by wavelengths in the millimetric infrared range, on top 
of which one observes rotational bands of energies growing as J(J -I- 1), with J non-negative integer, 
which disintegrate by emitting nanometric microwave radiations (so called rovibron model [T]). For 
a fixed vibrational mode, the microwave radiation is entirely attributed to the proportionality of the 
molecular Hamiltonian to the squared angular momentum operator, J^, 
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where /i stands for the reduced mass of the two molecules, and Re denotes the equilibrium bound length. 
The description of the rotational modes amounts to the diagonalizing of 'HrotiS, 7^) in the basis of the 
functions, $(d, (f) of the polar, 9, and azimuthal, </?, angles, 

'Hroti0,V’M0,P) = E^9,ip). ( 2 ) 

Interpreting as the moment of inertia, I, about the center of mass, the equation @ is cast in the 
standard form of the linear rigid rotor, termed to here as rigid rotator, according to, 

BJ^{9,7^)Yf{9,7,) = EjYfi9,ip), Ej=BJ{J+l), B = ^, (3) 

with B being the rotational constant. The solutions to the latter equation are the well known three- 
dimensional spherical harmonics, Yf^{9,ip), i.e. ^{9,(f) = Yf^{9,(p), defined as, 

Yf{9,7^)=E^‘^Pf{cos9), J = 0,1,2,... (4) 

Here, Pj^{cos 9) are the associated Legendre functions [5], known to become the Legendre polynomials, 
Pj(cosd), for M = 0. It is of common use to cast the rotational spectrum in Q in terms of frequencies, 
i/(J), as 

HJ) = ^Erot = l(Bj- = 2^J, (5) 

and define spectral lines separated by the constant gap of 2B/h. The above considerations refer to 
the idealized case of a pure rotator, the reality being that the equidistant line spacings are altered by 
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centrifugal distortions and other secondary effects, left aside in what follows. The rotational spectra 
are remarkable through the fact that the energy, Ej, equivalently, the frequency, in ([S]), depends 

on J alone meaning that the multiplicity of states in a level is (2J + l)-fold. Also the shapes of the 
disintegration (“radiation”) patterns, more precisely, the angular probability density distributions, given 
by the spherical polar plots of the squared {6, (p) functions, are energy degenerate. The plot in fig. 
[T] illustrates the degeneracy of the probability density distributions responsible for the disintegration of 
the states |J = 1,M = 0,±1 > to the ground state. 

Our principal goal is to answer the question whether the shapes of the probability density distributions 
are uniquely determined by the spectral frequencies, which parallels the celebrated question asked by 
Mark Kac in 1966, on whether the sound frequencies of a drum specify its form in an unique way [3]. 
Our case is that similarly to Kac’s question, also the answer to the one posed here by us turns out to be 
negative. 

Our point is that the perturbation of by a properly designed potential can retain the spec¬ 

tral frequencies and alter the wave functions. In this fashion, two sets of identical isospectral frequency 
patterns will correspond to two sets of probability density distributions of non-equivalent shapes. Stated 
differently, according to our findings, (2 J -|- l)-fold degeneracies do not necessarily imply equality of the 
Hamiltonian neither to the one of the canonical rigid rotator, nor to a similarity transformation of it. 
To be specific, below we show that the energies in the spectrum of the rigid rotator, perturbed by the 
trigonometric Scarf potential, remain identical to those within the spectrum of the free rotator, though 
the wave functions of the former are of unspecified parities and quite distinct from those of the latter, 
which are the spherical harmonics of well defined parities. The consequence will be rotational levels 
of non-vanishing electric dipole moments. The paper is structured as follows. In the next section we 
place the goal of our investigation within the context of spectral geometry and isospectral potentials. 
In section 3 we introduce the perturbation of the rigid rotator by the trigonometric Scarf potential and 
present its spectrum and wave functions. In section 4 we analyze the wave functions emerging after 
the perturbation and show they have lost their spatial parities, thus allowing for non-vanishing electric 
dipole expectation values. The paper closes with brief conclusions. 

2 Isospectral potentials from the perspective of spectral geometry 

Spectral geometry studies the relation between spectral properties of elliptic differential or pseudo¬ 
differentials operators and the associated Riemannian manifolds, among them curved surfaces. A Rie- 
mannian manifold, {M,g), is characterized by its metric g with covariant components, gap, and its 
Laplacian which is a second order differential operator acting on symmetric traceless tensors of sec¬ 
ond ranks, viewed as infinitesimal metric transformations. The general definition of a Laplacian in the 
spirit of Lichnerowicz reads, = —div(V) = — V“Vq = with the minus sign ensuring a 

spectrum bounded from below. 

A further key characteristic of a curved surface, in the following restricted to a compact Riemannian 
space, is the eigenvalue problem of the Laplace operator, Ag on this metric, 

Agfk = Xkfk, f = 0,1,2,.... (6) 

with fk being the corresponding eigenfunctions. The discrete number sequence, (0 =)Ao < Ai < A 2 < 
... — >■ 00 is referred to as the spectrum of Ag, or, of M. The vibrations of a thin membrane (a drum) are 
fixed by the eigenfunctions of the relevant Laplacian with Dirichlet boundary conditions and the square 
roots, '/Xi, of the eigenvalues are proportional to the sound frequencies. The fundamental questions 
asked by mathematicians in spectral geometries (see [4] for a pedagogical text) concern the properties 
of the spectrum given the geometry (the metric) and vice verse, whether the geometry (the metric) 
can be restored given the spectrum. While the first question is straightforwardly answered by solving 
the eigenvalue problem of the Laplacian, the answer to the inverse spectrum problem is far from being 
obvious. The issue is that the information provided by the spectrum alone is insufficient to allow one to 
uniquely recover the metric of the surface. Mark Kac’s question mentioned in the introduction on whether 
one can hear the shape of a drum [5] expresses precisely this very non-uniqueness of the solutions of the 
inverse spectral problem. Indeed, several drums of strictly equal spectra and non-equivalent shapes have 
been constructed beginning with work by Gordon and Wilson in 1984 [^. The applications of spectral 
geometry to physics are highly interesting and allow one to handle the so called isospectral potential 
problems [6]. The latter refer to two flat space Schrbdinger equations with distinct potentials, Vi and 
^ 2 , 
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which nonetheless describe exactly same spectra. Here, M denotes the mass of the particle under con¬ 
sideration. Several such potential pairs are known in quantum mechanics, the Morse and the hyperbolic 
Scarf potentials being one of them (see [7] for details). In one of the possibilities, though not necessarily, 
two isospectral potentials can be related by a so called isospectral transformation which is no more than 
a similarity transformation of any one of the Hamiltonians, say by some properly selected invertible 
smooth function, F, according to 
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meaning 
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(9) 


In general, potentials generated in this manner are complicated and contain gradient terms. The link 
of the Schrddinger equation to geometry (so called “geometrization” of Schrodinger’s equation 0 ) is 
established on the basis of the observation that the equations in © equivalently rewrite to eigenvalue 
problems of Laplacians in the so-called Maupertuis-Jacobi metrics, here denoted by (14) and defined 

as. 


g\f'\Vk) = = 2M{E, - Vk), fc = 1, 2, 


( 10 ) 


as 


^2^(m.)(^,) =H-2v4 (11) 

Here, is the so called “conformal” factor. Correspondingly, the metrics g\^'^\Vk) are referred to as 
conformal, while g(Mj )are the conformal Laplacians on the surfaces with the Maupertuis -Jacobi 
metrics, g^^'E{yfS) (with k = 1,2). It should be obvious that the equations in (ITT]) refer to isospectral 
problems on surfaces of distinct metrics and can be treated with the methods of spectral geometry, as 
indicated in the opening of this section. 

Back to the inverse spectral problem, in order to restore its uniqueness, new information needs to be 
invoked. As such, mathematicians consider the numbers of nodes and the nodal lines of the eigenfunctions 
[S]. Within the context of spectral geometry the goal of the present study can be cast in the following 
way: 


We draw attention to the fact that the angular potentials 
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( 12 ) 


are isospectral. 

However, rather than explicitly constructing the respective Maupertuis-Jacobi metrics, we here, being 
focused on physical applications in general, and on molecular physics in particular, prefer to visualize 
the difference of the aforementioned metrics indirectly. In first place the distinction is made through 
the evidently different shapes of the Vi{9) and 14(0) probability density distributions corresponding to 
equal quantum numbers and same energy, and in second place through the unspecified parities of the 
14(0) eigenfunctions, in contrast to the well defined parities of the 14(0) eigenfunctions. These are all 
physical entities of high spectroscopic relevance and are of major interest to experimental studies of the 
disintegration modes. 

Another interesting question regarding isospectral potentials concerns their symmetry properties. As 
we shall see in the following, the 14(0) potential is rotationally invariant because in properly chosen 
coordinates, it can be viewed as part of the Laplace operator on the two-dimensional spherical surface, 
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. The symmetry properties of the second potential are not instantly obvious and need to be elaborated 
separately. For this purpose one first has to systematically search for a complete set of operators which 
commute with j) and and then investigate their commutators with the aim to possibly figure 

out a relevant Lie algebra. This question, in going far beyond the scope of our investigation, will be left 
aside for the time being, and will be attended properly in a future work. The next section is devoted to 
the presentation of our case. 

3 Exotic degeneracy preserving perturbation of the rigid rotator 

In the current section we choose to simplify notations for the sake of transparency of the formulas and 
work in dimensionless units, setting B = 1. In eq. the parametrization of the two dimensional sphere 
is such that the polar angle 0 has been read off from the positive z axis (the North pole). Instead, we here 
take the liberty (without any loss of generality) to read off this angle from the negative y-axis, meaning 
that instead of tracing meridians in the usual way along North-East-South-West-North direction, we 
trace them along West-North-East-South-West. The reason behind this re-parametrization of the sphere 
will be explained shortly below, after the equation m- In so doing, the expression for and the 
spherical harmonics correspondingly changes to. 
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(13) 

(14) 

(15) 


With that in mind, we now consider the perturbation of the rigid-rotator by the following potential, 

v(0), 


3\9,ip) ^ H{9,ip)=3\9,ip)+V{9), 

6^ 26Mtan0 1 
cos^ 9 cos 9 4 

The eigenvalue problem of the perturbed Hamiltonian now takes the form, 

H{9, ^)^{9, = evl/(d, 4/(0, ^) = 


(16) 


(17) 


The second order operator o contains a first order derivative, which can be eliminated under the 
following variable change, 

U{9) = ^f^(j){9). (18) 

In effect, m acquires the shape of the following one-dimensional (ID) Schrodinger equation in the 9 
variable. 


n{9)U{9) = eUiO), n{9) 
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with the Schrodinger perturbation potential, here denoted by Vsci(9), being given as. 


6^ -I- a(a -1-1) b (2a + 1) tan 9 1 

cos^ 9 cos 9 4 ’ 

\M\ - i. (20) 

The squared magnetic quantum number, M^, in the operator on the sphere in (1131) . gets replaced in the 
one-dimensional Schrodinger operator by, 

m 2 ^ i = a(a+ 1) = (^\M\ - (^|M| + , (21) 

and acquires meaning of an ordinary potential parameter. 

According to the nomenclature in [7] , the singular interaction, Vsd (9) in (12011 is known under the name of 
the trigonometric Scarf potential, abbreviated, “ Scarf I”. This nomenclature, though strictly speaking at 
variance with Scarf’s original work |10] , where a csc^ singular potential has been considered, pays tribute 


Vsci{9) = 
a = 
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to the sech^ interaction, commonly named in the literature as hyperbolic Scarf potential, to which it 
relates by a complexification of the argument directly and without the need of introducing an additional 
shift by 7r/2. Instead, the csc^ potential, according to same nomenclature, is referred to as Poschl-Teller 
potential with emphasis on m where a combined csc^ + sec^ interaction has been introduced for the first 
time. Admittedly, other nomenclatures exist |12j according to which csc^ is referred to as trigonometric 
Scarf potential. As it will become clear in due course shortly below, after eq. (I25L the nomenclature 
does not affect the results. 


The singular 'H{0) eigenvalue problem is known to be exactly solvable and has been well elaborated 
in the literature. The general method consists in finding a suited point-canonical transformation that 
reduces the Schrodinger equation to the hyper geometric differential equation, known to have three 
regular singular points m, namely at x = 0, x = 1, and at infinity. At regular points, the solutions 
can be expanded in Frobenius series which in the case under consideration happen to terminate and 
give rise to Jacobi polynomials. This method has been pioneered by Stevenson as early as 1941 in 
[l4] . Alternatively, such equations can also be handled by the so called factorization method underlying 
modern super-symmetric quantum mechanics (see [7] for details), or by employing the recently elaborated 
Nikiforov-Uvarov method m- 

The solutions to (HU) are given, among others, in m- Upon equipping the wave functions with 
quantum numbers according to, 

[/(0)^[/l“l(0), (22) 

they are expressed as, 

= W(6») cos'^l+i 6»pl,“l-*'’l^l+^(sin6»), 

b 

W{6) = (23) 

\l — Sin 6 J 

Ct = tit + 1), t = \M\ -\- n, \M\ G [0, f], t = 0,1, 2,..., (24) 


and with n standing for the degree of the Jacobi polynomial, P“’^(sin6*) [TB]. Setting 6 = 0 recovers 
(modulo the VcosO factor in (fT^ l the solution of the free rigid rotator by equivalently re-expressing the 
spherical harmonics in terms of the Jacobi polynomials according to. 


'/cosO 


Pt“(sin6») = cos'll 0pl,^l’l^l(sin6l), n = t- \M\. 


(25) 


The 6 = 0 case corresponds to a Schrodinger equation with the sec^ 9 potential Vi (9) in (TT^ , while 
6 7 ^ 0 corresponds to (EHl), same as V2{9) in (fT^ . As long as the energy in (l24l) is 6 independent, the 
isospectrality of the potentials in (HU) is revealed. Notice that because the polar angle in o has been 
read off from the negative y-axis instead as usual from the North Pole, the argument in the associated 
Legendre functions is cos (f — 6*) = sin0 instead of the usual cos 9. In reality, we are dealing with Jacobi 
Polynomials of an argument x varying in the interval, — 1 < x < 1, and both choices for x, x = sin0, 
with 9 G [—7r/2, -|-7r/2], or, x = — cos6> with 9 =G [0, tt] are equally valid. 

Accordingly, the parities of the wave functions, to be studied in the next section, are obtained from 
9^—9 and not from 9 ^ tt — 9. Finally, the M values are restricted to natural numbers, satisfying the 
conditions, 

|M|-6+l>0, |M|-b6+l>0, (26) 

as required to ensure the known orthogonality of the Jacobi polynomials on the interval — 1 < sin 0 < -fl. 
As far as the \M\ values are integer, this restriction implies |6| < 1 for \M\ = 0. 


Comparison of (l24)) to S-dl, with both t and J taking same non-negative integer values, reveals 
the indistinguishability of the energy spectra of the free rigid rotator and the one perturbed by the 
trigonometric Scarf potential in the equations, (fTCl) ~ (l^^ . However the Scarf I wave functions in (1^51) are 
distinct from the spherical harmonics in (1151) . as it should be, given the non-commutativity of the squared 
angular momentum operator with this very potential, i.e. [J^(0, </?), Vsc/(0)] ^ 0. There is indeed no 
theorem which prevents two non-commuting operators of having equal eigenvalues. The operators non¬ 
commutativity only prohibits their simultaneous diagonalizing in same function basis. 

Our findings on the exotic rotational bands parallel at the quantum level the non-equivalent 

isospectral classical drums [3] mentioned in the introduction. 
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The next section is devoted to the loss of parity of the wave functions of the perturbed rigid rotator 
under investigation. 

4 Parityless states and non-vanishing electric dipole moments 

In the present section we study the relationship between the wave functions of the free and perturbed 
rotators. We begin with the simplest case of the ground state, t = 0, implying, n = M = 0. The 
corresponding wave function is read off from as 

Uo{S) = a/cos 0W(0) = \/cos0(l + sin 0)5 (1 — sin0)“^. (27) 

As a technical detail, we wish to notice that the ground state wave function, ipQ{9,ip), acts as the 
eigenfunction to an so(3) algebra Casimir invariant in a representation non-equivalent to the canonical 
one and defined by the following similarity transformation of the canonical squared angular momentum 
operator, J^, in (IT^ : 


3^{9,^)Yf{9,^) 

Yfi9,ip) 


cosi 9W{9)J^{9, ip)W-\9) cos"i 9, 
JiJ + l)Yf{9,^), 
V^W{9)Yf{9,ip). 


Indeed, the 3'^(9,ip) ground state, Y§{9,‘p), satisfies 


3\9,^)YS{9,^) = Q, 


(28) 


(29) 


whose obvious solution, 


YS{9, ip) = cosi 0W(0)roO(0, (^), Fo°(0, p) 



(30) 


coincides (modulo multiplicative constants) with ifjQ{9,p) in (l?7l) . The eigenfunctions, Y^[9,p), 
are obtained by rescaling the ordinary spherical harmonics, Yj^{9,p), by the transformation function, 
\/cos 0W(0), and will be termed to as “rescaled harmonics”. However, due to the unspecified parity of 
the function, W(0) in (1^51) . which is transformed by reflection into its inverse, W(— 6*) —^ W“^(0), rather 
than into itself up to a sign, all the rescaled harmonics are of unspecified parity. Related examples can 
be found in m- 


Our next example concerns the first excited state, corresponding to t = 1, and to the maximal allowed 
\M\ value of \M\ = 1, i.e. n = 0. One finds, 

^1(9,p) = Ul{9)e^‘^ = W(6i)cos5 9p-^'+^[si-a9)e^^ = W(6i)cos2 9P^ 

= y/^'W{9)Pl{siiv9)P^ =Y^{9,p), (31) 

and again an eigenfunction to the so(3) Casimir invariant in (12811 . 

3\9,p)Y^{9,p)=2Y^{9,p). (32) 

In general, it can be shown that all wave functions of maximal \M\ = t values (they include the ground 
sate) behave as J^(0, p) eigenfunctions and are characterized by a single angular momentum value, same 
that defines their energy. For all these states one is allowed to maintain the t label as a number conserved 
under transformations of the so(3) algebra in the representation defined in (l2^ . despite the loss of the 
parity specification explained after the equation (I30|) above. 

Next we consider consider the case of t = 2 and \M\ = 1, which corresponds to n = 1. We again express 
the relevant wave function %l)\{9,p) in terms of associated Legendre functions, Pj^^ {sin 9). In so doing, 
one finds 


4’\{9,p) = y/ cos0W {9) cos^’^^^(sin6()e*‘^ 
= Vcos 9'W{9) cos 9 [—&-I-2 sin6(] 
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= Vcos9W{e) 
= Vcos0W(6») 

= vW0W(6») 


[—& cos 9 + 2 cos 9 sin 9] 
-6pi(sin6»)e*‘^ + ^P2^(sin0)e*‘^ 


= -hYl{9,^) + ‘^-Y^{9M- 


As long as the parity of the spherical harmonics is (—1)'^, the functions Y^{9,ip)^ and Y2{9,+) are of 
opposite parities, and ip2{9,+) is built on top of a parity-mixed state. Moreover, the mixture of the 
rescaled harmonics in (1551) no longer behaves as eigenstate, and the t = 2 label in this case can 
not be interpreted as a quantum number conserved under transformations of the so(3) algebra in the 
representation of (l29ll . Decompositions of the art presented in (1551) can be performed for all the higher 
t values. They conserve the \M\ label and mix (rescaled) spherical harmonics with J = \M\, 




^ cjYf{9,+). 


J=\M\ 


(34) 


As a reminder, the Y^{9, cp) functions, defined in (1551) . are always of unspecified parity. It is obvious that 
the electric dipole moment, whose selections rules are AM = 0, and AJ = I, will have a non-vanishing 
expectation value in the tp2{9,(p) state due to the allowed Y2{9,(p) —>• Yi{9,(p) transition, on the one 
side, and more generally, due to the parityless nature of the W(d) function, on the other side. We have 
checked that also the general integrals 

de = J [ip^ {9, (p)Y sin 911;^ {9, (fi) sin 9d9d(p 0, (35) 

are non-vanishing for b values satisfying the conditions in (l26l) . We recall that the argument of the 
associated Legendre functions in these spherical harmonics is sin6>, versus cos 9 in the standard ones. 
Finally a comment on the symmetry properties of the wave functions of the perturbed motion is in order. 
These can vary and are determined besides by the \M\ value, as explained in the paragraph after eq. (1551) . 
also by the values of the b parameter. While for continuous b values the wave functions with \M\ ^ t do 
not behave as eigenfunctions, for regular integer, or half-integer b values, call them b = M', it can be 
shown that the Jacobi polynomials times cosl^^l 9 would behave as Wigner’s small functions [18] 

and the rotational symmetry is respected. We here emphasize on continuous b values. Comparison of 
the ordinary spherical harmonics, the wave functions of the free rigid rotator, to the Schrodinger wave 
functions of the Scarf I- perturbed rotator, is presented in fig. [5| 

5 Conclusions 

In conclusion, pure rotational frequency patterns do not necessarily imply equality of the Hamiltonian 
to that of the canonical rigid rotator. 

For discrete values of the potential parameter b, such patterns occur because the isospectral potentials 
in HU) are related through a similarity transformation of the type given in ([5|)-([5|) with the place of 
F being taken by the cos 2 0W(6)) factor in the wave function in (1551) . We discussed this issue after 
eq. dsg). Stated differently, the Schrodinger Hamiltonian with the trigonometric Scarf potential for 
such parameter values is identical to the similarity transformed squared angular momentum operator, 

= COS 2 dW(6))J^ cos “2 9'W~^{9), as explained in (1551) . Same happens when the magnetic quantum 
number \M\ takes its maximal value. In all these cases the reason behind the rotational patterns of the 
perturbed motion is still the so(3) symmetry of the Hamiltonian, though in a representation distinct from 
the canonical one. However, for continuous b values the above considerations are no longer valid and 
the reason behind the rotational spectrum is not obvious. In order to get an insight into the symmetry 
properties of the perturbed Hamiltonian for this case one has to follow the standard procedure and find 
the complete set of operators with which it is commuting. In calculating the commutators of these 
operators, a clue on a possible Lie algebra might be obtained m- 

In all cases, the wave functions of the rigid rotator, perturbed by the trigonometric Scarf potential in 
(EOD, despite of the rotational multiplicity of the states in the levels, are in first place always of unspecified 
parities. In consequence, a sort of exotics is produced in so far as the states in the rotational bands acquire 
non-vanishing electric dipole moments. The symmetry behind the emerging exotic rotational spectrum 
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for continuous b values is not very clear so far. We conclude on the possibility for existence of exotic 
rotational bands containing parityless states in two-body systems of a possibly new class, be them di¬ 
atomic-, di-molecular, quarkish, etc. and whose states carry non-vanishing electric dipole expectation 
values. 
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(a) |yO(0,v5))P {h)\Yt\9,^)\^ 

Figure 1: The three different shapes of the probability density distributions according to which one and 
the same energy, like Ej=i = 2B, corresponding to the first excited level, is released in the |1M >—>• |00 > 
ground state transition, i.e. spherical polar plots of the 11"° P and {9,ip) p functions, defined 
in O- The right figure counts twice because the radiation can go clockwise or, counter-clockwise. 
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Figure 2: Illustration of how equal rotational energies (frequencies) can correspond to non-equivivalent 
shapes of the probability density distributions. For the canonical rigid rotor, a probability density 
distribution of the release to the ground state of the energy, Ej = Et = 6B, of the second excited level, 
is given in fig. 2(a), while that corresponding to the perturbed rotator is displayed in fig. 2(b). The 
squared canonical rigid rotator wave function in fig. 2(a) is single J- valued, while the perturbed one in 
fig. 2(b) contains a AJ = 1 mixture according to ([321) ■ The potential parameter b has been given the 
generic continuous value of, b = 0.45. The figure visualizes the statement that one can not ’’hear” the 
shape of the radiation patterns from the spectral frequencies. 


9 






